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The deformability of a compact object induced by a perturbing tidal field is encoded in the tidal 
Love numbers, which depend sensibly on the object’s internal structure. These numbers are known 
only for static, spherically-symmetric objects. As a first step to compute the tidal Love numbers of 
a spinning compact star, here we extend powerful perturbative techniques to compute the exterior 
geometry of a spinning object distorted by an axisymmetric tidal field to second order in the angular 
momentum. The spin of the object introduces couplings between electric and magnetic deformations 
and new classes of induced Love numbers emerge. For example, a spinning object immersed in a 
quadrupolar, electric tidal field can acquire some induced mass, spin, quadrupole, octupole and 
hexadecapole moments to second order in the spin. The deformations are encoded in a set of 
inhomogeneous differential equations which, remarkably, can be solved analytically in vacuum. We 
discuss certain subtleties in defining the tidal Love numbers in general relativity, which are due to 
the difficulty in separating the tidal field from the linear response of the object in the solution, even 
in the static case. By extending the standard procedure to identify the linear response in the static 
case, we prove analytically that the Love numbers of a Kerr black hole remain zero to second order 
in the spin. As a by-product, we provide the explicit form for a slowly-rotating, tidally-deformed 
Kerr black hole to quadratic order in the spin, and discuss its geodesic and geometrical properties. 

PACS numbers: 04.20.-q, 04.25.-g, 04.70.Bw, 04.30.-w. 


I. INTRODUCTION 
A. Context of this work 

The relativistic theory of the tidal deformations of a 
compact object and the dynamics of a tidally-distorted 
self-gravitating body are fascinating and challenging 
problems in general relativity, which have received con¬ 
siderable attention in recent years. In a binary system at 
large orbital separation, the tidal interaction is negligi¬ 
ble and the two objects can be treated as point particles. 
However, as the orbit shrinks due to gravitational-wave 
emission, tidal interactions become increasingly impor¬ 
tant, and deform the multipolar structure of each object, 
their gravitational field, and the orbital motion up to the 
merger, where tidal deformations are dramatic and can¬ 
not longer be treated as small perturbations. Such tidal 
effects can leave a detectable imprint in the grav¬ 

itational waveform emitted by a neutron-star binary in 
the late stages of its orbital evolution, this system being 
the main target of current-generation gravitational-wave 
detectors [l-@. 

The multipole moments of a compact object in a bi¬ 
nary system are deformed by the tidal field produced 
by its companion. In Newtonian gravity, the constants 
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of proportionality between the multipole moments of a 
mass distribution and the perturbing external tidal field 
in which the object is immersed are known as tidal Love 
numbers (see e.g. [B,Q)- These numbers depend sensibly 
on the object’s internal structure and can thus provide a 
mean to understand thephysics of neutron-star cores at 
ultranuclear density Motivated by the prospect of 

measuring the Love numbers through gravitational-wave 
detections and also by the need of improv¬ 

ing gravitational-wave templates (e.g. [2l|, cf. 
for a review), in recent years a relativistic theory of 
tidal Love numbers has been developed with consider¬ 
able success [2, [23 - 1^ . The tidal Love numbers of static 
neutron stars have been computed for various realistic 
equations of state with great precision , whereas 

Schwarzschild black holes (BHs) were shown to have pre¬ 
cisely zero tidal Love numbers and, therefore, the 
multipolar structure of a Schwarzschild BH is unaffected 
by a perturbing tidal field. Very recently, this result was 
extended beyond the perturbative level 

These investigations have experienced a second burst 
of activity since the discovery of a set of relations among 
the moment of inertia, the electric tidal Love number, 
and the spin-quadrupole moment of a neutron star, which 
were remarkably found to be almost independent of the 
equation of state [1^ H^. These relations have been 
studied and extended in various scenarios [Sll - I^ and 
allow us to compute two elements of the triad once a 
single one is measured. For example, a single measure¬ 
ment of the tidal Love number of a neutron star would 
allow us - provided independent measurements of the 
mass and spin are available - to compute the moment of 
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inertia and the spin-quadrupole moment. Further, uni¬ 
versal relations have also been discovered, relating the 
mass quadrupole moment, the current octupole moment, 
and higher-order multipole moments (4^144^ . 

A further motivation to study tidal interactions in 
gravitational systems comes from the relativistic dynam¬ 
ics of an extreme-mass ratio inspiral, i.e. a stellar-size 
compact object being captured by a supermassive BH, 
which is one of the main targets of future space-based 
detectors [l^ . In the late stages of the inspiral before the 
plunge, the small object can perform hundreds of thou¬ 
sands of cycles around the supermassive BH, deforming 
the geometry of the latter by an amount proportional to 
the mass ratio This translates into changes to 

the orbital motion and to the gravitational-wave phase 
which are comparable to self-force [Hi effects. 

So far, the relativistic tidal Love numbers have been 
considered only for nonspinning objects. This is clearly 
a strong limitation, because compact objects in binary 
systems are expected to possess a nonvanishing angular 
momentum, and also because the coupling between the 
spin and the tidal field can produce novel effects that 
have been neglected in previous works. 

In this paper, we start a long-term effort to compute 
the tidal Love numbers of a spinning compact object and 
to study the tidal deformability in the presence of spin. 
Some complementary aspects of this complex problem 
have been investigated in the past. The geometry of a 
tidally-deformed, spinning BH was obtained in Ref. m 
by integrating the Teukolsky equations, whereas Ref. 541 
studied the intrinsic geometry of a spinning BH distorted 
by a small compact companion in the extreme-mass ratio 
limit. The latter two papers considered arbitrary values 
of the BH spin parameter. Very recently, Poisson com¬ 
puted the geometry of a slowly-spinning BH distorted by 
a quadrupolar electric and magnetic tidal field to first 
order in the BH spin [s^. He decomposed the tidal 
perturbations in terms of irreducible potentials and ex¬ 
tracted the tidal moments through a matching to the 
post-Newtonian metric of a binary system. Although 
perturbative in the spin, the work of Ref. [55j | provides 
physical insight on the spin-tidal coupling and presents 
the deformed solution in a simple and elegant way. 


B. Executive summary 

Here we summarize our main results and compare our 
work with previous analyses. Given the complexity and 
richness of the problem, the purpose of the present pa¬ 
per is manifold. First, we wish to develop a framework 
to compute the tidal Love numbers of spinning compact 
stars and we present here the first necessary step in that 
direction. This requires going beyond the analyses of 
Refs. [HI - IHj |. which only considered tidal deformations 
of spinning BHs. Similarly to Ref. (55j |. our framework is 
that of relativistic perturbations of a spinning geometry 
in the slow-rotation approximation [developed in Sec. [IT] 


and in Appendix |B] . Our results are valid to second or¬ 
der in the object’s angular momentum. A perturbative 
expansion in the spin is a reasonable assumption, first 
because neutron stars in binary systems are expected to 
rotate slowly, and also because the quadrupole moment 
(which describes the most relevant tidal deformations) is 
quadratic in the angular momentum. 

We find the general solution of Einstein’s equations de¬ 
scribing - in a closed analytical form - the exterior of a 
slowly-spinning object deformed by a tidal field generated 
by distant sources. Both the inner object and the distant 
sources are stationary and axisymmetric, and Einstein’s 
equations are linearized in the perturbations induced by 
the tidal field. Since the general solution is a linear com¬ 
bination of independent solutions, it depends on a num¬ 
ber of arbitrary constants. The values of these constants 
should be fixed by matching with the solution describing 
the interior of the compact star, or by imposing regu¬ 
larity conditions at the horizon if the central object is a 
BH. 

The structure of the general solution before fixing the 
constants allows us to discriminate the terms of the 
spacetime metric describing the tidal field from those de¬ 
scribing the linear response of the compact object to the 
latter. This discrimination is crucial to properly define 
the mass and current multipole moments [^, (M^ 

and respectively) of the central object in the buffer 
zone p, [2J, . Once the multipole moments are de¬ 

termined, the Love number can be obtained as the coeffi¬ 
cients of proportionality between the multipole moments 
and the tidal field, i.e. we define the quantities 


. (M) _ dMf, (5) _ dSi 


( 1 ) 


where £q is the amplitude of the tidal field. Because 
we include only linear corrections in the tidal field, the 
quantities above are independent of ^q. 

As we discuss in Sec. IHID H discriminating between 
the tidal field and the linear response of the compact ob¬ 
ject is a delicate issue, especially in the rotating case. 
Indeed, there is always the freedom to incorporate part 
of the response solution into the external tidal field, thus 
shifting the values of the multipole moments and then 
of the Love numbers (see e.g. [H, [HI - IHI| 1. To fix this 
ambiguity, here we follow a prescription, which we be¬ 
lieve is reasonable and it is also consistent with previous 
literature on the subject [1, HI, HU. However, this 
prescription is more well founded when the inner object 
is a BH than when it is a star. A more complete inves¬ 
tigation of this delicate issue is left for a future publica¬ 
tion; in particular, it is important to understand whether 
the Love numbers obtained with this choice correspond 
to actual measurable quantities, such for instance those 
appearing in post-Newtonian waveforms. 

For our general solution, we extract the Geroch-Hansen 
multipole moments (Sbl . [57j | and compute the formal ex¬ 
pressions for the tidally-induced changes of the first rel¬ 
evant mass and current moments from Eq. m- As dis- 
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cussed in detail in Sec. im the spin of the object intro¬ 
duces couplings between electric and magnetic deforma¬ 
tions and new classes of induced Love numbers emerge 
[cf. Eqs. (I5(l|) - (l55]) ]. In a future publication, the results 
presented here will be used to explicitly compute the tidal 
Love numbers (|5nil - (l55|) of a spinning neutron star to sec¬ 
ond order in the spin. 

When the object is a BH, the metric simplifies con¬ 
siderably. In this case, we impose the boundary condi¬ 
tions by requiring that the metric is regular at the hori¬ 
zon. The explicit solution that we present in Sec. m 
and Appendix |F] extends that derived in Ref. to 
include second-order effects in the spin but - at vari¬ 
ance with the solution of Ref. - is limited to the 
axisymmetric case and to the case in which the magnetic 
quadrupolar component of the tidal tensor is zero. As 
a by-product of our analysis, we compute a novel family 
of tidal Love numbers of the Kerr BH. Using the sep¬ 
aration of tidal and response solutions discussed above 
[cf. Sec. IIIID 1| . we find that the multipole moments 
of a tidally-distorted Kerr geometry are unaffected by 
the tidal source to second order in the spin and there¬ 
fore the corresponding Love numbers (HD vanish [see dis¬ 
cussion in Sec. IIV A| . This result extends the work of 
Ref. [ 25 , which found that the tidal Love numbers of a 
Schwarzschild BH are zero [bSlpl . Working independently 
from us, Landry and Poisson have recently obtained a 
complementary result [^ . namely they found that the 
Love numbers of a Kerr BH are zero to first order in the 
spin but, differently from ours, their results also include 
nonaxisymmetric tidal perturbations and a quadrupolar 
magnetic tidal field. 

Finally, our analytical BH solution can be used to com¬ 
pute easily the near-horizon deviations from a Kerr ge¬ 
ometry induced by the external tidal field. In Sec. IlYl we 
provide analytical expressions for the epicyclic frequen¬ 
cies, curvature invariants and other geometrical quanti¬ 
ties which can be relevant to investigate the dynamics of 
a binary system containing a spinning BH. We conclude 
in Sec. |Y|by discussing various interesting extensions of 
our work. 


C. Working assumptions 

The relativistic dynamics of a spinning object im¬ 
mersed in a tidal field is a challenging problem and we 
make a number of working assumptions to treat it. For 
the reader’s convenience, we list these assumptions be¬ 
low: 

(i) The central object is slowly spinning, with dimen¬ 
sionless spin parameter x = J/M^ <C 1, where M 


^ Very recently, the no-hair properties of tidally-distorted 
Schwarzschild BHs have been proved also beyond the pertur¬ 
bative regime [SSl . 


and J are the object’s mass and angular momen¬ 
tum, respectively. Our analysis is valid to 0{x^). 
Neglecting rotation and tidal deformations, the cen¬ 
tral is object is spherically symmetric. 

(ii) The tidal field varies slowly in time, so that time 
derivatives in Finstein’s equations are small com¬ 
pared to spatial derivatives. In practice, we consider 
stationary tidal perturbations. 

(iii) The object’s spin axis coincides with the axis of 
symmetry of the tidal field, which is assumed to be 
axisymmetric. This prevents the tidal field and the 
spin axis from processing [^, thus leaving the ge¬ 
ometry stationary. By adopting this assumption, 
here we mostly focus on the axisymmetric case, 
m = 0, where m is the azimuthal number of the tidal 
perturbations. Our method might be extended to 
the nonaxisymmetric case (describing, for instance, 
the complete tidal held generated by an orbiting 
companion), although precession introduces time 
dependence in the problem. We note that, because 
the background is axisymmetric, perturbations with 
different values of m are decoupled from each other. 
Therefore, our approach can describe the axisym¬ 
metric tidal response of the object, regardeless of 
the symmetries of the external tidal held. 

(iv) The sources of the tidal held are localized at large 
distance tq. Our vacuum solution is valid in the 
region i? < r <C cq, where R is the radius of the 
object. 

(v) The tidal held is weak: our results are valid to linear 
order in the amplitude of the tidal held. Assuming 
the tidal held is generated by a ring of mass me [as 
to satisfy the hypothesis of axisymmetry of point 
(iii) above], our results are valid for mcM'^/r^ <C 1. 

(vi) For simplicity, we assume that the tidal held is 
mostly electric and quadrupolar and consider 1 = 2 
polar tidal perturbations in the static case, where 
I is the standard harmonic index. This assumption 
is not crucial and can be easily lifted by including 
different t multipoles and a magnetic component. 

Through this work, we use G = c = 1 units. 

II. FRAMEWORK: SLOW-ROTATION 
EXPANSION 

The linearized dynamics of a spinning perturbed object 
in general relativity is a challenging problem because of 
mode mixing in the perturbation equations. However, a 
perturbative expansion in the angular velocity of the ob¬ 
ject can render the problem tractable. Our framework is 
that of linear perturbations of a slowly-rotating star, that 
has been initiated in Refs. (hSl - l^ and recently extended 
and put on hrmer basis in the context of BH perturba¬ 
tions 00 (see [t^ for a review). 
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The basic idea is that slowly-rotating geometries are 
“close enough” to spherical symmetry that an approxi¬ 
mate separation of the perturbation equations in radial 
and angular parts becomes possible. The perturbation 
functions are expanded in spherical harmonics and they 
reduce to a system of differential equations where various 
couplings between different multipolar indices (. and be¬ 
tween perturbations with different parity are introduced. 
The slow-rotation approximation imposes some selection 
rules on couplings between different multipoles and only 
a certain number of higher multipoles contributes to a 
given order in = Q/Hk <C 1, where is the uni¬ 
form angular velocity of the star and = \/M/R^ is 
the mass-shedding frequenc}0. This makes the method 
well suited to investigate complicated systems of coupled 
equations, which can be finally integrated by standard 
methods 

We consider a stationary, external tidal gravitational 
field which is characterized in terms of the Weyl’s ten¬ 
sor Cijki evaluated in the local rest frame, £ij = Ciojo 
and Bij = ^(ipqCjQ for the electric and magnetic tidal 
quadrupole moments, respectively. Because of our work¬ 
ing assumptions, it suffices to consider only stationary 
and axisymmetric perturbations of the spinning back¬ 
ground, although our technique can be generalized to dy¬ 
namical perturbations of any spinning compact objects 
to arbitrary order in the spin. 

A. Spinning background 

Following Hartle and Thorne [ 73 , the most gen¬ 
eral stationary and axisymmetricQ metric to 0{e^) 
in rotation can be written as 

ds'^= —e'^ [1 -I- 2e^ (jo + 72 ^ 2 )] dt^ 

, 1 + 2e^(mo -b m^Pi)! {r - 2M) ^ ^ 
l-2Mlr 

+r‘^ [1 -b 2el{v2 - j2)P2] [dd^ + sin^ - CaUjdtf] , 

( 2 ) 

where P 2 = P 2 (cos'd) = (Scos^d— l)/2 is a Legendre 
polynomial. The radial functions v and A4 are of ze¬ 
roth order in rotation, w is of first order, and jo, j 2 , mo, 
m 2 , V 2 are of second order. We also introduce the func¬ 
tions e~^ = 1 — 2A4 /r and uj = fl — uj at zeroth and 
first order in the rotation, respectively. By plugging this 


^ When considering BHs, the expansion parameter is naturally 
ea = X- Hereafter we use Ea as a bookkeeping parameter for 
the expansion in the angular momentum. 

® We also require the spacetime to be symmetric with respect to 
the equatorial plane, and invariant under the “circularity condi¬ 
tion”, t —>■ —f and ip —>■ —ip, which implies gt^ = gtq, = gr-B = 
3r-(p = 0 I75I . Note that, while the circularity condition follows 
from Einstein and Maxwell equations in electrovacuum, it might 
not hold true for other matter fields. 


decomposition into the gravitational equations = 0 
and by solving the equations order by order in e^, we ob¬ 
tain a system of ordinary differential equations (ODEs) 
for the rotating background [zl, 13 . To second order 
in the spii 0 , the vacuum background metric is given in 
Appendix This solution describes the exterior of a 
spinning object, with mass M in the static case and di¬ 
mensionless spin parameter x <C 1. As discussed below, 
the constants O, 5m and 5q which appear in the solution 
are related to the angular velocity of the central object 
and to 0 (x^)- corrections to the mass and and to the 
quadrupole moment, respectively. 

The geometry given in Appendix lAl describes a slowly- 
rotating Kerr BH as a particular case. In this case, 
regularity of the metric at the Schwarzschild horizon 
r = 2M imposes Sq = 0, whereas 5m can be set to 
zero without loss of generality through a mass rescaling. 
Our coordinates differ from the more standard Boyer- 
Lindquist ones. In the metric ([2|) the horizon’s location 
is r+ = 2M [1 — e^x^/ 8 ], whereas the ergoregion is lo¬ 
cated at Tergo = 2M [1 — e^x^ cos(2i?)/8]. Although the 
coordinates are singular at r = 2 M, all curvature invari¬ 
ants are regular outside the horizon. 

B. Perturbations 

Slowly rotating and oscillating compact objects can be 
studied as perturbations of the axisymmetric, station¬ 
ary solutions discussed above. For completeness, here we 
discuss the case of nonaxisymmetric perturbations with 
azimuthal number m and harmonic index £ > \m\. Later, 
we shall restrict to m = 0. Scalar, vector and tensor field 
equations in the background metric ([ 2 |) can be linearized 
in the field perturbations. Any perturbation function 5X 
can be expanded in a complete basis of spherical harmon¬ 
ics; schematically, we have 

, (3) 

where we have imposed that the perturbations do not de¬ 
pend on t explicitly and is a basis of scalar, vector 

or tensor harmonics, depending on the tensorial nature 
of the perturbation 5X . As in the spherically symmetric 
case, the perturbation variables 5Xg^ can be classified 
as “polar” or “axial” depending on their behavior under 
parity transformations. 

The linear response of the system is fully characterized 
by a coupled system of ODEs in the perturbation func¬ 
tions 5Xg^^. In the case of a spherically symmetric back¬ 
ground, perturbations with different values of (£, to), as 
well as perturbations with opposite parity, are decoupled. 
In a rotating, axially symmetric background, perturba¬ 
tions with different values of to are still decoupled but 
perturbations with different values of i are not. 


This expansion has been recently extended to 0{e‘^) l45l . 
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To second order, the perturbation equations read 
schematically (see Ref. [t^ for a pedagogical derivation) 

0 = A-i + eaTuAe + e^Ai 


+ ia{Q{Pt-l + Qn+lPt+l) 

+ Cq Q(,-iQiAt-2 + Qi+2Qi+iAt+2 +C>(eDi(4) 


Since in the nonrotating limit axial and polar perturba¬ 
tions are decoupled, a possible consistent set of solutions 
of the system (HI)-® has = 0 = 0, 

where .^ = L is a specific value of the harmonic in¬ 
dex. This ansatz leads to the “axial-led” [t^ subset of 
Eqs. 


^ = 1^1 + eamVi + AgVe. 

+ ^a{QiAi-i + Qe+iAe+i) 

+ e-a Qi-iQi'Pe-2 + Qi+2Qi+i'Pe+2 

where we have defined 
Qe = 


£2 _ jj£2 

4£2 - 1 ’ 


0{el), (5) 


( 6 ) 


and At-, An, An, An, At are linear combinations of the 
axial perturbations with multipolar index t, similarly, 
'Pe, Pe, Pe, Pe, Pe are linear combinations of the polar 
perturbations with index i. 

The structure of Eqs. is interesting. In the 

limit of slow rotation a Laporte-like “selection rule” 
imposes perturbations with a given parity and index i 
to couple only to: (i) perturbations with opposite parity 
and index £ ± 1 at 0{ea)', (h) perturbations with same 
parity and same index ^ up to 0{e^)\ (hi) perturbations 
with same parity and index £± 2 at C*(e^). Furthermore, 
from Eq. ® it follows that Q±m = 0, and therefore if 
\m\ = £ the coupling of perturbations with index £ to 
perturbations with indices £ — 1 and £ — 2 is suppressed. 
This general property is known as “propensity rule” 
in atomic theory, and states that transitions £ —>■ £-1-1 are 
strongly favored over transitions £—>■£— 1. Note that 
the slow-rotation technique is well known in quantum 
mechanics and the coefficients Qt are in fact related to 
the usual Clebsch-Gordan coefficients [t^ . 


C. Axial-led and polar-led perturbations 

Due to the coupling between different multipolar in¬ 
dices, Eqs. 0-® form an infinite system of coupled 
ODEs and the spectrum of their solutions is extremely 
rich. However, in special configurations the perturbation 
equations can be greatly simplified, as we now show. 

First, we expand the axial and polar perturbation func¬ 
tions (schematically denoted as atm and ptm, respec¬ 
tively) that appear in Eqs. dl]) and ®: 

atm = + Ca 42 + + C’(ea) 

Ptm =42 + Ea 42 + 442+^>(4) • (7) 

The terms At ±2 and Pi ±2 in Eqs. 0-® are multiplied 
by factors e^, so they only depend on the zeroth-order 

perturbation functions, afl^m^ 4±2m- The terms At±i 
and Vt±i are multiplied by factors ta, so they only de¬ 
pend on zeroth- and first-order perturbation functions 
( 0 ) ( 0 ) ( 1 ) ( 1 ) 

Pt±lm’ ^t±lm’ Pt±lm- 


Al + eatnAL + ^aAr + f^a^QbPL-1 + Ql+iPl + i) = 0 
'Pl + 1 + ia'mpL + 1 + iaQb+lAL = 0 

< 'Pl-1 + f-a^npL-l + eaQhAh = 0 

Al+ 2 + iaQL+2pL+l + £0 2^+12i+2“4L = 0 
. Al-2 + ^aQL-lpL-1 + f^aO-LQL-lAh = 0 

( 8 ) 

where the first equation is solved to second order in the 
spin, the second and the third equations do not contain 
zeroth-order quantities in the spin, and the last two equa¬ 
tions do not contain zeroth- and first-order terms in the 
spin, i.e. aL ±2 = The truncation above is consis¬ 

tent because in the axial equations for £ = L the polar 
source terms with £ = L± 1 appear multiplied by a factor 
Ca, SO terms Pj^^im would be of higher order in the axial 
equations. 

Similarly, another consistent set of solutions of the 

same system has p 2 ± 2 m = ^ “iiim = The corre¬ 
sponding “polar-led” system reads 

'Pl + f^a'mpL + + f-aiQbAL-l + Ql+iAl + i) = 0 

Al+1 + eatnAL+i + eaQh+iPh = 0 

< Al-1 + ea^AL-1 + eaQLpL = 0 

Pl+ 2 + ^aQL+2AL + l + PiQL+lQL+2pL = 0 
. 'Pl-2 + CaQi-i-Ti-i -I- CqQlQl-i’Pl = 0 

(9) 

Interestingly, within this perturbative scheme a notion of 
“conserved quantum number” L is still meaningful: even 
though, for any given L, rotation couples terms with op¬ 
posite parity and different multipolar index, the subsys¬ 
tems ® and ® are closed, i.e. they contain a finite 
number of equations which fully describe the dynamics 
to second order in the spin. 

The main assumption that leads to Eq. ® [resp. 
Eq. ®] is that only axial (resp. polar) perturbations 
with harmonic index L are activated at zeroth order in 
the rotation. In terms of an external tidal field, we are 
assuming that such field is a pure £ = L magnetic (resp. 
electric) state at zeroth order in the rotation. This as¬ 
sumption would not hold if the zeroth-order tidal held 
is a mixture between different £ states. In such case 
one has to deal with the full system ®-®, which is 
much more involved. However, working with the sys¬ 
tem ® should provide a reliable approximation, because 
the electric quadrupolar (£ = 2 ) contribution to the ex¬ 
ternal tidal held is the dominant one. 

The explicit form of the axial-led and polar-led sys¬ 
tems ® and ® for a spinning stationary and axisym- 
metric object is derived in Appendix [B] and is available 
in a Mathematica® notebook provided in the Supple¬ 
mental Material. 
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III. TIDAL DEFORMATIONS TO 
SECOND-ORDER IN THE SPIN 

Let us now focus on axisymmetric (m = 0) stationary 
perturbations, which are a particular case of those dis¬ 
cussed above and in Appendix |B] The axial-led (resp. 
polar-led) system ([8]) [resp. dH)] describes the second- 
order spin corrections to a nonrotating object immersed 
in an £ = L magnetic (resp. electric) tidal field, which 
we now compute. For clarity, we focus on the dominant 
L = 2 polar-led perturbations, although the same proce¬ 
dure can be applied to other values of L and to axial-led 
perturbations. As discussed above, L = 2 polar pertur¬ 
bations in the static case couple to L = 1 and L = 3 axial 
perturbations to first order in the spin, and to L = 0, 
L = 2 and L = 4 polar perturbations to second order 
in the spin. Within this perturbative scheme, the per¬ 
turbation equations at each order are naturally written 
as inhomogeneous ODEs, where the homogeneous part 
depends on differential operators defined in the nonro¬ 
tating case, whereas the sources depend on the couplings 
between different perturbations. For this reason, it is use¬ 
ful to present the equations in the static case for generic 
values of L. 

Unless otherwise written, henceforth we reabsorb the 
bookkeeping parameter Cq in x- 


A. Tidal perturbations at zeroth order 


In the static case, axial and polar perturbations as well 
as perturbations with different harmonic indices are de¬ 
coupled from each other. To zeroth order in the spin, the 
axial and polar sectors of stationary perturbations with 
i = L > 2 reduce to two single decoupled homogeneous 
ODEs @,[23: 




r(r-2M)“° 

L{L + l)r(r - 2M) + ,l) 


VAAhf^]=hf^" + 


r2(r — 2MY 
AM - rL{L + 1) 


HA’ = 0 ( 10 ) 


r^ir — 2M) 


= 0, 


( 11 ) 


where we have defined the differential operators Dp,l and 
'Da,l- The metric function is algebraically related 
to and its derivatives, whereas the functions 
and vanish at zeroth order in the spin and therefore 
the polar perturbations are also static. Focusing for sim¬ 
plicity on the L = 2 polar sector, the explicit solution 


reads @,[13 


^(2) _ 2>ar{r - 2M) ^ 7 

0 - m2 +2 


3r(r - 2M) 


M2 


log 1 


1 


2M / 6 

-h r I- 

r \M r — 2M 


- 5 


: 


= 3a 


M2 


- 2 


/ 3r^ 

I 2M2 


- 3 log 1 


3r 

'M 


2M 


-b3 



( 12 ) 

(13) 

A) 

(14) 


where a and 7 are two integration constants, so that the 
solution above is a linear combination of two independent 
solutions. 

The solution proportional to a diverges at large dis¬ 
tances and can be identified with the external tidal 
field @,113, HI, whereas the solution proportional to 7 
naturally represents the object’s linear response to the 
applied tidal perturbation (see however Sec. IIII D II for 
a more detailed discussion). As we shall discuss later, 
the constant 7 is associated with the tidally-induced 
quadrupole moment of the object in the static case, 
whereas a is proportional to the axisymmetric com¬ 
ponent of the electric quadrupolar tidal field (roughly 
speaking, a ^ rricM^/r^ for a ring of mass rric at orbital 
distance rg). The ratio y/g is proportional to the L = 2 
electric Love number @, |24| , 


(2) _ A-fM^ 
15ai?5 ’ 


(15) 


where R is the radius of the object. The constants a and 
7 are proportional to each other, their ratio being de¬ 
termined by matching the exterior solution above to the 
regular solution describing the perturbed object’s inte¬ 
rior. Thus, for an object of given mass and composition, 
( 2 ) 

the Love number ' is uniquely determined. In the 
static case, the dependence of the Love number on the 
stellar equation of state has been discussed in detail in 
Refs. @,@,i[l3,[23,[2i,p. 

In case the central object is a BH, regularity of the 
geometry at the Schwarzschild horizon r_|_ = 2M im¬ 
poses 7 = 0, as can be directly seen from Eq. (1121) or 
by computing some curvature invariant at the horizon. 
Therefore, the Love number (1151) of a tidally-deformed 
Schwarzschild BH is precisely zero. This is a general re¬ 
sult, which is valid for any L and for electric and magnetic 
tidal perturbations [23. In other words, the multipolar 
structure of a static BH is not deformed by a perturbing 
tidal field (very recently, this result has been extended to 
arbitrary values of a static tidal field [^1. 


B. First-order corrections 

The zeroth-order solution (HHl-dll sources the axial 
perturbations with L = 1 and L = 3 through the second 


















7 


and third equations in the system ([9]), yielding 


VA,2[hl^^] = 4'^" 


4(3r-M) ^( 3 ) ^ ( 3 ) 
r2(r-2M) ° ^ 


(16) 

(17) 


where 'DA,i are differential operator^ and are source 
terms that are given in Appendix [Cj The axial met¬ 
ric functions and vanish identically. As ex¬ 
pected, the sources are proportional to the coupling be¬ 
tween the background gyromagnetic term, gt,p, and the 

zeroth-order function . It is easy to verify that the 
first-order corrections to the polar perturbations are van¬ 
ishing, so Eq. m and the two equations above fully 
characterize the polar-led L = 2 system to first order in 
the spin. The explicit solution of the equations above is 
given in Appendix [Dl This solution depends on four new 
integration constants, 01,3 and 71,3, which arise from the 
homogeneous problem associated with Eqs. (flBll and (fT71) . 
and that are discussed in Sec. mmi 


C. Second-order corrections 


With the zeroth-order and the first-order solutions at 
hand, from the first equation in the system ([9]) we can 
compute the second-order correction to the metric coeffi¬ 
cient which we denote by 6HQ^\r) to distinguish 

it from the zeroth order quantity. This correction satis¬ 
fies the following inhomogeneous ODE 


Vp^2[6H^''^] = 


2{r-M) (2)' 

r{r-2M) ° 


2{2AP + 3r^-6Mr) 


r2(r — 2M)2 


= Sp 


( 2 ) 


(18) 


where the source Sp'^ is also given in Appendix [Cl Note 
that SHl ' = 0 [i.e. polar perturbations remain static 
also to 0{x'^)], whereas and are algebraically 

related to and its derivatives. 

Finally, to fully characterize the second-order correc¬ 
tions, one needs to compute the last two equations in ©, 
which define the second-order terms in the induced L = 0 
and A = 4 polar sectors. The A = 4 system reduces to 
the second-order ODE 




2 (r-M) ( 4 )^ 

+ r(r- 2M)‘^ 0 


4 (m 2 -P 5r2 - lOA/r) ( 4 ) _ 


r2(r-2M)2 


= Sp 


(4) 


(19) 


^ Note that, while 'Dj\ 2 is obtained from the operator ^ for 
L = 2, the operator i is different because L = 1 pertur¬ 
bations satisfy a different set of equations, as discussed in Ap¬ 
pendix 


where the source S'p ^ is given in Appendix |C] Also in 
this case the other A = 4 polar components follow alge¬ 
braically from 5H^'^ and its derivatives. 

On the other hand, the A = 0 polar system satisfies a 
different set of equations (cf. Appendix |B|), which can be 
reduced to the following first-order system: 




r-2M 
^ r-2M 


c.(o,o) 
Op , 

c.(0.2) 

Op , 


( 20 ) 

( 21 ) 


and the sources Sp’^'^ and Sp’^^ are given in Appendix (Cl 
Remarkably, all the equations above can be solved an¬ 
alytically. Schematically, the nonvanishing metric coef¬ 
ficients to quadratic order in the spin read (reinstating 
the bookkeeping parameter Ca only in these equations) 


gtt = -e" 


l + 2el{jo + j2P2-^-^{n-Cof 


+ (fff + ) 

gtip = - w) sin^ d 

+e,sin^(4'yi° + /i("y3o) , 


Y 


20 


Qrr — 


1 - 


2M 


' „ 2 Wq -f TO2P2 

r-2M 


+elSH^^^Y^ 


t{0)^00 


( 


H. 


( 2 ) 




gm = r2 [1 -p 2 e^(u 2 - j 2 )P 2 
+el5K^°'^Y°° + ( a :(2) 


2MA(2)) 


Y 


-20 




) p. 


( 22 ) 

(23) 


(24) 


(25) 

(26) 


where we recall that = y^°(d) are the scalar spher¬ 
ical harmonics with m = 0 and P 2 = 2^/tFJSY‘^^ is a 
Legendre polynomial. The radial functions v, Ai, lu, jo, 
J2 j 477.0 j 477.2, V 2 are given in Appendix [^ the radial func¬ 
tions and are given in Appendix |D1 the radial 
functions Hq^\ and Ar^2) are given in Eqs. (IT^ - ([T4)l : 
whereas the radial functions SHq^\ 5H^\ 

5P[^\ are cumbersome and, 

to avoid typographical errors and help comparison, their 
full solution is provided in an online notebook in the Sup¬ 
plemental Material. Note that the only nonvanishing off- 
diagonal term of the metric is gt,p and it only contains 
the background gyromagnetic term and the axial pertur¬ 
bations with A = 1 and A = 3. 





















D. Description of the solution 


Equations (l2^ - (l26l) fully describe the exterior metric 
of a tidally-deformed spinning object to second order in 
the spin and in the region R < r ^ ro. Because the pro¬ 
cedure to obtain such solution is considerably involved, 
as a nontrivial consistency check we have verified that the 
explicit solution satisfies Einstein’s equations in vacuum, 
Rfii' = 0, to quadratic order in the spin and to linear 
order in the tidal perturbations. 

The background Hartle-Thorne solution depends on 

I 


the parameters M, y, fl, 6m and Sq; the 0(x°) tidal so¬ 
lution depends on the constants a and 7; the 0(x) tidal 
solution depends on the constants ai_3 and 71,3; finally, 
the 0{x^) tidal solution depends on the constants 00,2,4 
and 7o,2,4- The full solution depends on 17 free param¬ 
eters, whose physical meaning is summarized in Table H] 
and will be discussed in more detail below. Note that 
each pair of constants and 7^ arises from the homoge¬ 
neous problem associated with the inhomogeneous equa¬ 
tions for the corresponding multipole i presented above. 

As a representative example, we present here the ex¬ 
plicit form of the function , which is the most com¬ 
pact among the second-order perturbations: 


= aox^ + 


8aiX , 27ix^ , 2 

+ + ax 


2 - y y/i{y - 2) V3(y - 2)y^ 


4(6 -f y{2y - 1)(1 -f y{6y - 1))) 
V5{y - 2)y^ 


( 3 (-4 -f y(-4 +{y- 2)y(3(y - l)y - 14)) -t 3(y - 2fy‘^{2 + y) coth ^ [1 - y]) 


+7X 


V^{y - 2)y 

36 - y(26 + 3y(2 + {y - l)y{15y - 11))) + §2/(24 + 2/(-12 + y{12 + y{-32 + (41 - IS//)//)))) log 


y-2 

y 


Sq 


2^5(2/-2) V 


3^(// - 2)//5 

(24 + 2y (-20 + y (256 + 9// (-14 + (// - 2fy))) 


+6(// - 2)2/(-4 -t 2/(-6 +{y- 2)y{3>{y - l)y - 14 ))) log 

[2/-2I 

+ ^(2/-2)V(2 + 2/)log 

[2/-21 
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[ y \ 

2 

L y \ 

)\ 


(27) 


where y = r/M. The parameters oq and 70 arise from 
the homogeneous problem associated with Eqs. (EOl)- 
m- The constant Oq only appears in gtt and, as we 
discuss below, can be eliminated through a time rescal¬ 
ing. The constants ai and 71 are related to the source 
terms proportional to the axial term with L = 1, hQ^\ 
which sources in Eqs. (I^ - (I^ through a coupling 
with the spin of the object, whereas the constants a and 
7 are related to the polar terms with T = 2 at zeroth 
order in the spin (namely and which 

couples to 6Rq^^ at second order in the spin. Clearly, 
the structure of the solution reflects the selection rules 
discussed in Sec. m 

More generically, the tidal corrections of the metric 
can be schematically written as a linear combination of 
independent solutions in the following form: 


= aSglZ^+xSgjZ^+Y^ [aeSglZ"'' + , (28) 


e=o 


where we have factored out the dependence on a, 7, ai 
and Xi- The first two functions Sgj^ and Sgj^J contain 
terms of zeroth order in the spin and also terms of first 
and second order in the spin which arise as the particular 
solutions of the inhomogeneous equations; the functions 


6g\Z‘\ 6gl2f'^ with odd i contain 0{x) terms arising from 
the solutions of the homogeneous problem and also O(x^) 
terms arising from the particular solutions of the inho¬ 
mogeneous equations; finally, the functions 5g\w \ 
with even £ contain only O(x^) terms. 


i. Separating the tidal part from the linear response 


Crucially, the only parts of the metric which diverge 
in the far-held limit are Sg^^ and SglZ‘\ whereas the 
terms Sg^^) and yield an asymptotically-hat solu¬ 

tion. However, note that one cannot simply identify the 
terms proportional to a and with those associated 
with the external tidal held, and those proportional to 
7 and Xi with the linear response of the central objects. 
Such dehnition suffers from an ambiguity because a triv¬ 
ial shift of the integration constants 


7 = 7 - 0^7 : 


a'^i 


(29) 
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TABLE I. List of the free parameters appearing in the tidally- 
deformed metric of a spinning vacuum geometry for polar-led 
L = 2 perturbations. The subscript in ae, and 7^ refers to the 
multipole that is related to the specific constant. The precise 
relation between the multipole moments and the constants 7 
and 7f is given in Eqs. (Slli-illHl). 


0{x,a) 


(0,0) 

M 

mass 

(1,0) 

X 

spin 

(1,0) 

Q. 

angular velocity of the object 

(2,0) 

5m 

spin-induced mass shift 

(2,0) 

5q 

spin-induced quadrupole-moment shift 

(0.1) 

a 

external electric quadrupolar tidal field 

(0,1) 

7 

static response to the external tidal field 

(1,1) 

ai 

external magnetic dipolar tidal field 

(1,1) 

71 

tidally-induced spin shift 

(2.1) 

OLQ 

constant shift of gtt at infinity 

(2,1) 

70 

tidally-induced mass shift 

(2,1) 

CX.2 

external electric quadrupolar tidal field 

(2,1) 

72 

tidally-induced quadrupole-moment shift 

(2,1) 

Q3 

external magnetic octupolar tidal field 

(2,1) 

73 

tidally-induced octupole-moment shift 

(2,1) 

Ct4 

external electric L = 4 tidal field 

(2,1) 

74 

tidally-induced hexadecapole-moment shift 


would transform Eq. (j28l) to the equivalent fornU 






e=o 


+ 1 ^ 912 ) + ■ 


(30) 


£=0 


While the new solution proportional to a is still divergent 
at large distances - and therefore a can still be identified 
with the amplitude of the tidal field - the coefficients 
of the subleading solutions and Sg^Z^ have been 

shifted by an amount proportional to a. In Sec. IIII El 
below we show that such coefficients are related to the 
multipole moments of the central object (but see discus¬ 
sion in Sec. IIIIFjl . It is therefore crucial to find a unique 
prescription to characterize the linear response of the sys¬ 
tem. 

In the nonrotating case, this ambiguity was mentioned 
in Refs. (UliSill- A rigorous way to identify the two 
solutions uniquely is to perform an analytical continua¬ 
tion in the number of spacetime dimensions d or in 
the multipolar index I [7^ and to recognize that the 


® As discussed below, a quadrupole-led electric tidal field at infinity 
imposes = 0. For clarity, we have already used this condition 
in Eq. (13011 and postpone its proof to Sec. IIII D 21 
^ We are indebted to Jan Steinhoff for suggesting Refs. [^ . |77H to 
us. 


general solution for the gravitational potential at large 
distances schematically reads [Hil 

(55«-^arqi + ...)+7r-^-‘'+3(l + ...) , (31) 


where the dots represent a series in Mjr which can also 
contain logarithmic terms and a and 7 are integration 
constants. Comparison with the Newtonian potential 
generated by an Cpole distribution allows us to identify 
the first and the second solution as those describing the 
tidal field and the linear response of the system, respec¬ 
tively. By treating d or £ as real parameters, the two so¬ 
lutions above can always be distinguished without mixing 
of the possible common powers in the series expansion. 
Such procedure give s a precise proof of the identification 
done in Refs. [2, [2^ for the nonrotating case. Indeed, it 
shows that the a- and 7-solutions appearing in Eq. dm 
precisely describe the tidal field and the linear response 
of the object, as anticipated. 

Unfortunately, performing a similar analytical contin¬ 
uation in the rotating case is impractical because the 
solution of Eq. m for generic L is written in terms 
of Legendre functions, which appear in the sources of 
Eqs. (HU)-(El). While, for any specific value of L, the 
latter equations admit a solution in closed form, the same 
does not seem true for generic L. Furthermore, the prob¬ 
lem gets only more involved to second order in the spin. 
For this reason, to separate the tidal part from the linear 
response in the general solution (|28|) we have adopted a 
different procedure. We fix the shifts in Eq. (|29|) such 
that the new growing solution after the shift. 


SgfZr^^ = a 



(32) 


only contains a finite numbers of terms in a large-distance 
expansion. As it turns out, when the central object is 
a BH this prescription selects a unique solution for the 
shifts (1^ , at least up to quadratic order in the spin. We 
claim that such solution represents the physical solution 
of the tidal field, whereas the remaining part represents 
the physical response of the BH to the tidal field. 

In the case in which the central object is a compact 
star, this truncation only occurs up to first order in the 
spin. Indeed, to quadratic order there appear infinite 
terms proportional to ax^Sq in the far-held expansion 
of which cannot be canceled by the shifts (1^ . 

This reflects a certain arbitrariness in the dehnition of the 
Love numbers of a spinning star, which we will investigate 
in the future. In this case, we hx the shifts (E^ such that 
the far-held expansion of the growing solution contains 
a hnite number of terms, modulo those proportional to 
ax^Sq, which arise from the particular solution of the 
inhomogeneous problem at second order in the spin. 

Although the one just described is admittedly not a 
rigorous prescription, it is nonetheless supported by the 
following observations: (i) it allows us to identify cor¬ 
rectly the two pieces of Eq. ED in the nonrotating case, 
in which the unique solution is known. Indeed, it is clear 
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from Eq. (IT^ that 7 = 0 is the only possibility, (ii) To 
first order in the spin, it agrees with the recent results by 
Landry and Poisson [63 that were obtained using an in¬ 
dependent prescriptioro. In particular, our prescription 
automatically implies that the tidal solution is regular 
in the exterior spacetime, also when the central object 
is a BH. This is due to the fact that possible terms in 
the form ^ log(r — 2M), ^ l/{r — 2M) that might ap¬ 
pear in are precisely removed using the prescription 
above. Finally (iii), our procedure incorporates the idea 
that we “superimpose” a well-behaved tidal field on a 
central object. The solution corresponding to the tidal 
field, regardless of whether the central object is a BH or 
a star, should not have a pole at r = 2M. We only allow 
poles arising from the coupling with irregular background 
terms, which are those proportional to ax^Sq discussed 
above. 

In the rest of this paper, we will identify the tidal solu¬ 
tion and the linear-response solution through the proce¬ 
dure just described and postpone a more rigorous analy¬ 
sis for future work. For simplicity, we rename the rescaled 
constants as 7 ' —?> 7 and 7 ^ —?> 7 ^. Note that the physical 
meaning of the constants listed in Table U refers to the 
solution after this rescaling. 

Finally, we observe that although the tidal solution 
extracted through this procedure contains only a finite 
number of terms in the series expansion in the BH case, 
it can nonetheless contain terms with the same powers of 
M/r as the linear-response solution. This is a peculiarity 
of the spinning solution and does not occur in the static 
case. For example, contains up to M/r terms 

(which also appear in 5g^/Z^) to first order in the spin, 
and up to (M/r)^ terms (which also appear in SgjZ^^) to 
second order in the spin. 


Boundary conditions at infinity 


The integration constants ai that appear in the vac¬ 
uum solution (1281) have to be fixed by some physical re¬ 
quirement on the nature of the tidal field. Because we 
assume an (L = 2 )-leading external tidal field, we fix the 
constants af, by imposing that the other leading L ^ 2- 
contributions vanish. 

Let us start by the first-order corrections in the spin. 
The large-distance behavior of the gyromagnetic term 
reads 


gtip -t asxM 


I 


5 


I 


y;^°sinz? 





(33) 


® We are indebted to Phil Landry and Eric Poisson for sharing 
their results with us which have helped in revising the argument 
above. 


The leading-order terms proportional to and 
would respectively correspond to some spurious L = 3 
and L = I magnetic tidal perturbations [ 2 ^ and we 
eliminate them by fixing as = 0 = ai. Note that the 
constants 71 and 73 do not appear in the leading-order 
behavior at large distance and are therefore not associ¬ 
ated to possible components of the tidal field. 

Likewise, to second order in the spin, the dominant 
large-distance behavior of the gtt metric component reads 

gtt ^ aiX^Y^°W^ + O . 


This leading term is related to a possible L = 4 compo¬ 
nent of the external electric tidal field, which we eliminate 
by fixing ^4 = 0. Again, the constant 74 does not affect 
this leading-order behavior and it is indeed related only 
to nondivergent terms. 

In a similar way, we fix a 2 by requiring that the 
leading-order behavior of the component gtt —t r^ at 
large distances is not affected by the spin. This has to be 
the case because the leading-order behavior is related to 
the components of the tidal field, which does not depend 
on the properties of the central object, as discussed be¬ 
low. Finally, the constant ao appearing in Eq. (|27| can 
be eliminated through a time rescaling t > (I -|- 
where ry is a constant to be fixed. To second order in the 
spin, this rescaling only affects the coefficient gtt and not 
the gyromagnetic term gt^p, because the latter is of the 
order %, so corrections would be cubic in x- Without loss 
of generality, we use this gauge freedom to fix 


2-yng = —ao + 



Sq, 


(34) 


which completely cancels the term ao in the metric and 
also simplifies the angular dependence of some subleading 
term in the gtt coefficient at large distances. Henceforth, 
we will enforce the time rescaling just described and that 
ai = 0 {£ = I,..,4). Thus, the free parameters of the 
tidal perturbations reduce to a, 7 and 7 ^ {i = 0, ..,4). 

To check the asymptotic behavior of the solution, let 
us identify the tidal field in terms of the metric compo¬ 
nents. First, we consider the large-distance expansion of 
our solution after hxing the constants a^. 


1 + ffti 3aV {'&) 


^ _ 4 (I + Smx^) ^ 


+4 (1 -I- 26mx^)] 


gtip -t -3\/ -ax sin f} 

\l ^ 


M 


3 -I- 5cos(2z?) 


grr — 1 —t 3aF + • ■ • ) 


^ - I ^ 3ar20(t?) 


^ -8v^(l-f 25mx^) 


-f ... 


together with the exact relation gpp = sin^ '&gm- In the 
expansion above, we neglected terms of the order Mjr or 
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higher, which depend on a, 7 and 7^. Note that the time 
rescaling discussed above has been used both to eliminate 
ao and to make the subleading, 0(r°), tidal term of gtt 
proportional to With a different gauge choice, 

the angular dependence would not factor out. 

From the large-distance expansion above we can iden¬ 
tify the tidal field. In a suitable gauge, the large-distance 
behavior of a vacuum spacetime distorted by a quadrupo- 
lar tidal field reads [1, [H, [bj 


gtt -1 - Sijx'x^ + ... , 

(35) 

gti ^ ^ 

(36) 

gij ^ ^ij [1 SijX 

(37) 

where is the Levi-Givita symbol, £ij and Bij are 

the electric and magnetic tidal quadrupole moments, re- 


spectively, and x* are Cartesian coordinates. We can 
now transform to spherical coordinates, x* = rn*, with 
n® = (siniJcos^, sinz?sin(/j,cosd), and decompose the 
tidal tensor as 13] 


2 

^ (38) 

m——2 

where are symmetric traceless tensors re¬ 
lated to the usual spherical harmonics by = 

By plugging this decomposition into 
Eqs. dsa-dsii), transforming the metric to spherical co¬ 
ordinates, and comparing with the leading-order asymp¬ 
totic behavior of our solution, it is straightforward to 
identify 

Sof 

£0 = ’ £m^o = 0, Bij = 0. (39) 

The gt^ coefficient in the large-distance expansion above 
contains only subleading terms which are not related to 
a magnetic quadrupole tidal moment, but are due to the 
frame dragging effect. The fact that the only nonvan¬ 
ishing component of the tidal field is £0 is a consistency 
check of our solution, since we imposed a pure electric 
and axisymmetric quadrupolar tidal field. For a tidal 
source of mass rric at a distance xq, £0 ~ rricIrQ [H, 
and Eq. (l3^ gives a ^ WcM^/xq, as previously antici¬ 
pated. 


E. Multipole moments 

In order to extract the multipole moments 0,113 from 
our solution, we have to remove the tidal fields, as in 
the nonrotating case n m m, 113. This is done by 
identifying the solutions corresponding to the tidal field 
and to the linear response of the system as discussed in 
Sec. IIII D II and by setting to zero the constants corre¬ 
sponding to the unique growing solution thus defined, 
i.e., by fixing a = ae = 0. As discussed in Sec. IIII D 11 


the remaining solution is asymptotically flat. In this wot, 
we determine the spacetime’s response to the tidal fielqj. 

We follow Ryan’s approach to compute the Geroch- 
Hansen multipole moments [s^, in a gauge-invariant 
way through the geodesic properties of an axisymmetric, 
asymptotically-flat spacetime [tqI - I^ . Ryan found that 
a low-velocity expansion of the energy change per log¬ 
arithmic interval of the orbital frequency is completely 
determined by the multipole moments of a Ricci-flat so¬ 
lution. Such quantity is defined as 

A£.-n^, (40) 


where G is the angular velocity of an equatorial circu¬ 
lar geodesic with specific e nerg y Ep. The large-distance 

m 


expansion of AE reads [4 
v'^ 20 Si 

y ~ y 

'225 Bos'll 70 M; 

ly ~ 




Mi 


28 Si 


27 Mi 

f 81 Si S1M2 

( „ ,.,+6 


9 Mi 

S. 


\ 2 Mi 

935 M 2 


-6 


Mi 


35 35 M4 


24 Mi 


12 Mi 


12 Mi 


y- 


,10 


/6615 115 Sf 


128 


18 Mi 


165 


Si 1408 968 S'! M2 


Mi 


243 Mi 


+ 


/ 45927 , 123 Sf 
256 
SiM2 


27 
9147 M2 


Mi 


^2 S 3 

~9 


Mi 


,,11 


14 Mi 


56 


Mi 


93 M| 

TMf 


-24 


MJ 




where v = (Moll)^A is the linear velocity, Me are the 
Geroch-Hansen mass multipole moments, whereas Se are 
the Geroch-Hansen current multipole moments [H, [13 . 
Note that these moments are equivalent [13 to the mul¬ 
tipole moments defined by Thorne [13 using asymptoti¬ 
cally mass-centered Gartesian coordinates, and that the 
definitions above are all gauge invariant. 

For the asymptotically-flat, stationary and axisym¬ 
metric spacetime discussed here and written as ds^ = 
gp^dx^dx'^, we have 


G = 

Ep = 


T y {gttp,r')‘^ gtt,rgtp(p,r 

9ip^p,r 

Qtt 

yj 9tt 


(42) 

(43) 


^ We remark that one cannot simply read off the multipole mo¬ 
ments from the components of the full solution, following e.g. 
the prescription of Indeed, this prescription requires that 

the metric is expressed in asymptotically mass-centered Carte¬ 
sian coordinates and, most importantly, that the spacetime is 
asymptotically flat. None of these conditions is satisfied by the 
full metric 1I28II . 
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By using the definitions above and inverting the func¬ 
tion V = v{r), we obtain the expression of AE given 
in Appendix [E) By comparing Eq. (lEll) with Eq. gB), 
we can identify the nonvanishing multipole moments to 
quadratic order in the spin: 


M 

A 

M 2 

M3 


A 

M4 

M4 

M5 


1 -f 6mx — 


2 7 o + 12 \/Zx 5 q 2 


1 -b 


4v^ 
y37i-2V57\ 
4v^ J' 


j 6 q-l]x^- 


27 


^/5^T 280^/^T - 

-b4V5(2872 -b 7(135 -b 56,5 to -b 765g)) 
3-\/773 -b 44\/57 


195^773 


-X, 


28^/tt 

65^773 -b 4 ( 42 O 74 - 757(221 -b 4325g)) 
294O0F 


-X 


(44) 

(45) 

(46) 

(47) 
2 

7 

(48) 


Note that, because the expansion (HTl) contains more co¬ 
efficients than multipole moments, the comparison with 
Eq. (HII) is also a nontrivial consistency check of our so¬ 
lution, because some relations between the coefficients of 
the small—u expansion are fixed through Einstein’s vac¬ 
uum equations. 

Interestingly - due to the coupling to L = 3 axial and 
L = 4 polar terms discussed in Sec. ITT] - the tidal field 
introduces a nonvanishing current octupole S 3 and a non¬ 
vanishing mass hexadecapole M 4 even if the background 
solution does not possess such moments to second order 
in the spin. Clearly, such corrections would add to the 
(spin-induced) S 3 and M 4 terms appearing in the back¬ 
ground solution to third and to fourth order in the spin, 
respectively (^ . However, since we neglect terms higher 
than second order in the spin, S3 and M 4 are absent from 
our background solution, and only appear in the tidally- 
induced deformation. This is similar to the nonspinning 
case, in which the tidal field induces a quadrupole mo¬ 
ment [the first term in Eq. (1461) ] even if the central object 
was originally spherically symmetric. 

The result above provides a practical definition of 
the tidally-induced multipole moments of the spacetime. 
Once the exterior metric is matched to the interior so¬ 
lution, the constants M, y, Sm, 6 q, 7 and 7 ^ can be 
extracted and can be related to the multipole moments 
of the external spacetime using the definitions above. 


We wish here to comment on some further technical is¬ 
sue which is related to the definition of the multipole 
moments. 

Our procedure is based on the fact that the perturba¬ 
tion equations are linear and, therefore, the two solutions 
mentioned above independently solve Einstein’s vacuum 
field equations. However, even if the full perturbed so¬ 
lution behaves - by definition - linearly, the multipole 
moments of the spacetime might be mixed among the 
two solutions. In other words, the multipole moments 
of the central object might in principle be contaminated 
by the external solution. An example of this fact is the 
static and axisymmetric Weyl solution, whose line ele¬ 
ment reads 

ds 2 ^ ^ ( 49 ^ 

where C/, k and W depend only on p and 2 . It is easy 
to prove that Einstein’s equations impose the potential 
?7 to be a harmonic function in flat space, V'^U = 0. 
Therefore, in this case a linear combination of two solu¬ 
tions (say Ui and U 2 ) is still a solution of the Laplace 
equation. However, the resulting moments will in gen¬ 
eral be a nonlinear combination of the moments of the 
individual solutions [1^. To linear order, Ui and U 2 can 
be considered as small perturbations and therefore non¬ 
linear terms can be neglected. However, this example 
shows that a mixing between two independent solutions 
can occur in the computation of the multipole moments. 

Eor this reason, it is not clear whether the moments of 
the linear-response solution are the true moments of the 
deformed compact object. We stress that this limitation 
is not a prerogative of our approach. The same criti¬ 
cism equally applies to the static case (see Ref. [ 2 ^ for 
a discussion). Although our procedure is reasonable and 
fully equivalent to previous approaches, we believe that a 
more rigorous analysis is needed to solve this important 
issue, even in the static case. This would likely require 
a fifth order post-Newtonian expansion of the field equa¬ 
tions for a binary system, which is currently not avail¬ 
able for generic mass ratios, in order to determine the 
quantities which actually appear in the post-Newtonian 
gravitational waveforms. Indeed, the definition of the 
tidal Love numbers in a relativistic theory would remain 
mostly academic without a proper connection to observ¬ 
able quantities. 

Leaving these problems for future work, in the follow¬ 
ing we simply follow the standard procedure to define the 
tidal Love numbers. 


F. Note on the definition of multipole moments G. Tidal Love numbers of a spinning object 


In Sec. IHIDll we have discussed some subtleties and 
some degree of arbitrariness in separating the solution 
describing the external tidal field from that describing 
the response of the system, even at the linearized level. 


It is clear from Eqs. (I44l) -(l48 p that the external tidal 
field modifies various multipole moments of the space- 
time. Because such corrections are necessarily linear in 
the tidal field £q, we define various tidal Love numbers 
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and as follows (cf. Eq. ([T|)): 


> (M) _ i^to m -r V uL/to I 2 

^0 -■ ’ 

,(S) V3d£,ji-2V5de„i 


9£o7o + l^V^dSfffSq 


(50) 

(51) 


.(M) _ 

A 2 — 


2 d£ffM^ 


l95V7d£,X3 


y/^ 280y/TT - 

+4'/5(289£(,72 + i9£o7(135 + 56i5to + 76dq)) 


(S) _ 3^a£o73 + 44y55f:o7 ^^4 


A 3 — 

AM) ^ 

2940v^ L 


280 F 

&7>V7d£Ai3 


M\, 


(52) 

(53) 

(54) 


-4 (420a5„74 + v^a£o7(221 + 432d9))] , (55) 


where Sfj, denotes derivative with respect to fo oc a [cf. 
Eq. (j39l) ]. Because the quantities above are linear in the 
tidal field, the derivatives with respect to £q are just num¬ 
bers (see Sec. IIII Hi belowi. The first term in Eq. (1521) is 
proportional to the standard electric quadrupolar Love 
number (1151) and it is the only term that does not de¬ 
pend on the spin. The other terms in Eq. (1521) , as well as 
Eqs. dsni)-® and Eqs. (ESI)-®, are novel spin-induced 
corrections. 

The tidally-induced corrections to the multipole mo¬ 
ments are linear- (or higher-) order quantities in the 
spin. For an electric quadrupolar tidal field, the mass 
quadrupole gets tidal-induced corrections at quadratic 
order in the spin, i.e. the corrections enter at the same 
order of the spin-induced quadrupole moment of the Kerr 
metric. On the other hand, both the current octupole S '3 
and the mass hexadecapole M 4 get tidally-induced cor¬ 
rections which enter at lower order in the spin than the 
spin-induced corrections, which would enter at 0 {x^) and 
0 {x‘^), respectively [isj] . 

More precisely, to quadratic order in the spin, the tidal 
correction to M 2 is suppressed by a factor a <C 1 rela¬ 
tive to the undeformed spin-induced term. However, the 
tidal corrections to and M 4 are dominant with re¬ 
spect to higher-order spin-induced corrections [ 1 ^ when¬ 
ever x^ a. This condition is consistent with our per¬ 
turbative scheme because both a and x ^^e small per¬ 
turbation parameters, although it might be difficult to 
match in practice, except for extremely slow rotations 
which would however make the spin-induced tidal Love 
numbers (I50D - E5P almost zero. 

Finally, here we focused on the most interesting case 
of the tidal Love numbers associated with a quadrupo¬ 
lar electric tidal field. Nonetheless, there is no reason 
to expect qualitatively different results for other compo¬ 
nents of the tidal field. The selection rules discussed in 
Sec.ini suggest that an axisymmetric electric tidal held 
with multipole (. would introduce corrections to the mass 
multipole moment Mi to 0 (x°) and O(x^), to the mass 


multipole moments Mi ±2 to O(x^), and would introduce 
corrections to the current multipole moments 5'^±i to 
0(x)- Likewise, we expect that a magnetic tidal held 
with multipole I would introduce corrections to the cur¬ 
rent multipole moment Si to 0 {x°) and 0 (x^), would 
modify the moments Si ±2 to O(x^), and would intro¬ 
duce corrections to the mass multipole moments Mi±i 
to 0 {x)- A natural extension of our work is to compute 
these novel families of Love numberJ^. 

In particular, an axisymmetric magnetic tidal held 
with £ = 3 should modify the mass quadrupole moment 
M 2 by a term linear in the spin and linear in the intensity 
of the tidal held. Because such term enters at linear order 
in the spin, it might be the dominant deformation of the 
mass quadrupole M 2 , although it would be suppressed by 
the fact that the octupolar magnetic component of the 
tidal held is much smaller than the quadrupolar electric 
component for typical sources. 


H. Extracting the Love numbers of a spinning 
object 

In general, computing the Love numbers (I5(I|) - (I55]) re¬ 
quires a numerical integration of the perturbation equa¬ 
tions in the interior of the central object and a matching 
procedure with the analytical exterior solution discussed 
here. 

If the central object is a self-gravitating huid, the in¬ 
terior solution consists of the metric perturbations and 
the huid perturbations, the latter vanishing in the ex¬ 
terior. By requiring regularity at the center of the ob¬ 
ject and continuity at the surface, the Love number (1151) 
in the static case can be extracted 0 from the ra¬ 
tio evaluated at the radius R of the star 

and by using the analytical solution m- Likewise, 
the constants (^ = 0,1,..,4) can be extracted from 

h^o\R) and SH^q'^\r), re¬ 
spectively, and using the analytical solution presented in 
this work. After this matching, the parameters 7 and 
Xi will be necessarily proportional to a oc fp and, there¬ 
fore, the dependence on £q of the Love numbers (l50ll - 
(1551) disappears, as expected. One can conveniently fac¬ 
torize the fg-dependence in the constants, by dehning 
say 7 o = and the matching procedure would 

then allow us to extract dmtidai, which represents a mass 
shift induced by the tidal field at quadratic order in the 
spin, similarly to the parameter 5m that represents a 
mass shift purely induced by the angular momentum of 


Note tliat an electric (resp. magnetic) tidal field with odd (resp. 
even) values of I would generate perturbations which break the 
reflection symmetry of the background Consistently with the 
argument just presented, such components of the tidal field would 
induce multipoles such as M 3 , S 2 , etc.., which are identically 
vanishing in the case of a solution with refiection symmetry as 
the one discussed here. 
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the background solution. Similarly, the tidal Love num¬ 
bers ((5IIll - (l551) also depend on the coupling between the 
static Love number oc 7 /a and the quadrupole shift 
dq, which also appears at 0 {ax^)- 

In a subsequent work, we will apply this procedure to 
compute the tidal Love numbers (1501) -(1 55 [) for a spinning 
neutron star. In the next section, we will instead focus 
on the case in which the central object is a BH, which 
can be remarkably solved analytically. 

IV. TIDALLY-DISTORTED SPINNING BHS 

The solution discussed above becomes much more 
tractable in the BH case. The latter can be obtained 
by requiring regularity of the solution across and outside 
the event horizon. Because the Hartle-Thorne coordi¬ 
nates [13. Eli in which our ansatz ([ 0 ]) is written are singu¬ 
lar at r = 2 M, the regularity of the metric perturbations 
is slightly more subtle than in a regular set of coordi¬ 
nates as, e.g., that adopted in Ref. [^. Nonetheless, to 
ensure regularity one can compute some curvature invari¬ 
ant, such as the Kretschmann scalar RabcdR°‘^'^‘^ and the 
Pontryagin density *RR = RabcdR^^'^ef, and impose 

regularity at the singular point. 

For the background solution given in Appendix |Al reg¬ 
ularity imposes Sq = 0 whereas dm can be set to zero 
without loss of generality through a redefinition of the 
mass of the background solution. Furthermore, as pre¬ 
viously mentioned, regularity also imposes 7 = 0 in the 
solution (|I2D - (II4|) . It is also easy to show that the cur¬ 
vature invariants are divergent at r = 2M unless 

72 = 73 = 74 = 0 . (56) 

Due to the parity properties of our solution, regularity 
of the Kretschmann scalar fixes the constants 72 and 74 
related to the even-parity (polar) perturbations, whereas 
regularity of Pontryagin density fixes the constant 73 as¬ 
sociated with the odd-parity (axial) perturbations. As 
shown in Eqs. (l4^ and (1^ . the constants 70 and 71 are 
related to tidally-induced mass and spin shifts, respec¬ 
tively, and can be set to zero without loss of generality 
in the BH case, because they can be reabsorbed in the 
definition of mass and spin of the background Kerr met¬ 
ric. After imposing 70 = 7 i = 0 and the conditions (1561) . 
the Kretschmann scalar reads 

PI^RabcdR^^^^ = ^ + [1 + 3cos(2iI)] 

yO V TT 2/*^ 

72y^ 

+ -^ (2 + 2/ - 6?/^ - (8y^ + y - 6) cos(2zI)) 

27qi 

+y(783 -k y(45y - 3II)))) - 4(104 -k y(-I04 
+y(48 + 5y(79 + 5(y - I)y)))) cos(2,?) 

-5(40 -k y(y(y(65 -k y(3I -k 35y)) - 56) - 48)) cos(4zI)] , 


and is regular everywhere except at r = 0. A similar 
regular expression can be obtained for the Pontryagin 
density. Note that the Kretschmann scalar can also be 
decomposed in spherical harmonics as M'^RabcdR°'^‘^‘^ = 
/ 2 i(?')y^*°('!^), where /o. 2 , 4 (?') are radial functions. 
The exact regular metric of a tidally distorted spinning 
BH to second order in the spin is given in Appendix [Fj 

As a consistency check, we have verified that, to first 
order in the spin, our solution reduced to that found 
by Poisson in Ref. (55j | in the axisymmetric case with 
zero magnetic quadrupolar component of the tidal field. 
The transformation between our coordinate system and 
the null coordinates of Ref. [11] is given in Appendix |D1 
Eqs. (lD3])-(lD6l). 

A. Love numbers of a Kerr BH 

We can now compute the tidal Love numbers for a 
spinning BH. By imposing regularity through (5y = 7 = 0, 
using Eq. dMI), and setting Sm = 70 = 7i = 0, from 
Eqs. (IT41) -(l48 p we obtain 

dM2 dSs dM^ 

dSo d£o d£o ' ^ ’ 

together with the equations = 0, which just 

represent the freedom of rescaling the mass and spin of 
the background metric. Therefore, we immediately ob¬ 
tain that the Love numbers of a Kerr BH are zero as in the 
Schwarzschild case [13, m and as recently found to first 
order in the spin by Landry and Poisson [b^. We note 
that the separation of solutions discussed in Sec. IHID II 
is crucial to obtain such result. Any other prescription 
would account for a shift as in Eq. (1291) which would mod¬ 
ify the Love numbers. Therefore, the prescription given 
above is the unique one that yields zeroth Love numbers 
for a Kerr BH to second order in the spin. On the light 
of these results, it is also natural to conjecture that the 
Love numbers of a Kerr BH are zero to any order, at 
least in the axisymmetric case considered here. 

As a by-product of the BH uniqueness and no-hair the¬ 
orems [HjII^ (see also [ssl - l^ l. the multipole moments 
of any stationary BH in isolation can be written as [s^j . 

Me + iSe = {ixY . (58) 

All moments with £> 2 can be written in terms of 
Mq = M and Si = J through the above relation. 
Therefore, any independent measurement of three mul¬ 
tipole moments (e.g. the mass, the spin and the mass 
quadrupole M 2 ) is a null-test of the Kerr metric and, 
in turn, it might provide the first genuine strong-gravity 
confirmation of general relativity [s^ - l^ . Our results 
(together with those of Ref. M) show that the no-hair 
relations (I58p are robust not only to nonperturbative ef¬ 
fects in the tidal field [13] , but also to perturbative tidal- 
spin interactions. In other words, the relation (1581) holds 
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also for a slowly-rotating BH immersed in a weak tidal 
field and it can therefore be interpreted as a general¬ 
ization of the no-hair theorems for stationary, tidally- 
deformed spinning BHs. 


B. Properties of the solution 


The BH solution given in Appendix [F1 is an analytical 
(albeit perturbative) solution of Einstein’s vacuum equa¬ 
tions and, as such, it is interesting per se. It is therefore 
relevant to study this solution more in detail, e.g. by 
computing various geometrical and geodesic quantities 
related to this spacetime. Let us start by computing the 
intrinsic geometry at the horizon. The horizon location 
is defined as the largest root of g'"’' and reads 


r+ = 2M 




(59) 


Although the horizon location does not depend on the 
coordinates and (/ 3 , its intrinsic geometry is not spher¬ 
ical. To compute the intrinsic metric, we consider the 
spatial section dt = 0 oi the metric ([2]) at r = r+. In the 
slowly-rotating limit, we obtain 

dSt=const,r=r+ = 9^^ (»' = »'+, "d)dn'^ . (60) 

and therefore the intrinsic geometry is spherical only 
when evaluated at r = r+ does not depend on d. 
Using the solution (IF4|) . the Ricci curvature of the in¬ 
trinsic geometry reads 

M^Rintr = ^ - cos(2z?) -f -a [1 -b 3cos(2zI)] 

Z O 4 V TT 

-= [1 + 172 cos(2d) + 115 cos(4i?)] , (61) 

128V 57r 

and is constant only in the nonrotating and undeformed 
case. The intrinsic curvature is shown in Fig. [T] for dif¬ 
ferent values of x and a. 

It is also interesting to look at the geodesic structure of 
the spacetime. For a stationary and axisymmetric space- 
time, one can define the potential 


V{r) = [Ep{r, 7 r/ 2 ) + l] , (62) 

where U{r, d) = 5 “ — 21 g*’^ + Pg‘^‘^, Ep is the energy per 
unit mass of a point particle (given in Eq. (1431) 1. and I is 
the proper angular momentum. For a circular, equatorial 
orbit at r = Tc, Ep and I can be determined by imposing 
f^(u) = Vi'fc) = 0 , where the prime indicates differen¬ 
tiation with respect to r. The solution for Ep is given 
in Eq. m- The further condition V" = 0 yields the 
location of the innermost stable circular orbit (ISCO), 
f = nsco- 

Furthermore, by considering small perturbations of a 
circular, equatorial orbit, one finds the epicyclic frequen¬ 
cies Hr and governing small oscillations in the radial 



FIG. 1. The curvature of the intrinsic geometry of a tidally- 
deformed spinning BH for different values of the spin x S'Hd 
of the intensity of the tidal field a. We used extreme values 
of X 8 -nd a in order to magnify the effect of the deforma¬ 
tions. The coordinates {X, Y) are related to (r, d) through 
r = VA'2-by2, Y/X = tan-i d. 


and in the d direction, respectively. These read 


= 

nl = 


{gtt + ^gtcp)'^ d'^u 

2 grr dr'^ 

{gtt + ^gt(p)^ 

2 gss dd^ 


(u,|) , (63) 

(r„f). (64) 


Using Eqs. (ga), dna and dH, it is straightforward to 
show that, in our case, the explicit form of the orbital 
and epicyclic frequencies reads 


MQ-—^ ^ 3 - 2) g X 

y3/2 Y TT 

(3 + ^)( 3 t/- 2 )x^ 

42/11/2 


Sox 

2 y/^y^ 
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{M^sY = — -J- 

yi 


5 3{y-2Y{2y-l)a 6x 

TT 


4y3 
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9/2 


case, 


3 (74 - 120y + 30y^ + ax (-6 + 25y + 9y^) 
^\f^y'^l’^ 2y'^ 

3 (-120 + 1048y - 948?/2 + 57y3 + 346?/'^) ax^ 

^ 

By evaluating these frequencies at the ISCO, we obtain 


^ISCO 


oISCO 

“i9,0 


1 + 0.150x0.2 + 0.022 xL + 0.017r3-o3 
-0.004r3-(fxo.2 + 0.00008r3-o3x2 2 , (71) 

1 + 0.123x0.2 + O.OlSxL + O.OlOr^o^ 
-0.003r3-(fxo.2 - 0.0002r3-o3x2 2 , (72) 
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(69) 


whereas the ISCO radial epicyclic frequency van¬ 

ishes to second order in the spin, as in the undeformed 
Kerr case. 


C. Tidal versus spin effects 


where each denominator denotes the corresponding fre¬ 
quency evaluated at x = ci = 0 , xo .2 = x/( 0 - 2 ) and 
^30 = ro/{30M). As expected, tidal corrections are small 
because the tidal field is suppressed by the third power 
of the orbital distance. However, for tq rs 30M (i.e., 
a ~ 10 “®) and x ~ 0 - 2 , the linear tidal correction is 
comparable to the quadratic in spin correction and they 
both correct the static, undeformed result by a few per¬ 
cent. 

With the normalization of x and cq adopted above, the 
correction linear in the spin and in the tidal field - which 
arises from the spin-tidal coupling - is of the order of 
0.5%, whereas the correction quadratic in the spin and 
linear in the tidal field is much smaller, of the order of 
0.02%. Note, however, that our metric is valid in the 
limit r <C To, whereas the frequencies above are evaluated 
at the ISCO, nsco 6 M, which is only moderately 
smaller than tq ~ 30M. For this reason, our discussion 
is intended only as an order-of-magnitude estimate. 


Our results can be used to estimate the effects of spin 
and tidal deformations on the BH geometry. The analysis 
of this section will be mostly qualitative, a more detailed 
study will appear elsewhere. Let us assume that the tidal 
field is generated by a source of mass rric at a distance 
ro from the central BlS. In this case, £q ~ TOc/i'o and 
Eq. (15^ yields 


which is understood as an order of magnitude estimate. 
We also assume rUc ~ M so that a ~ {M/tqY. In this 


We recall that we are considering an axisymmetric tidal field. 
Therefore, strictly speaking, our source should be a ring of mass 
rric and radius rg. Nonetheless, the qualitative results of this 
section would also apply to the more realistic case in which the 
source is a companion star of mass me at orbital distance ro. In 
this case / 0 S'Ud the tidal field sources both axisymmet¬ 

ric and nonaxisymmetric deformations in the metric. By virtue 
of the axisymmetry of the background, modes with different az¬ 
imuthal number m are decoupled from each other and the ax¬ 
isymmetric components of the metric perturbations are exactly 
described by our solution. 


V. CONCLUSIONS AND EXTENSIONS 

Computing the tidal Love numbers of a spinning neu¬ 
tron star is an open problem with various potential appli¬ 
cations in classical general relativity and in gravitational- 
wave astronomy. Here, we presented a framework to 
study gravitational perturbations of a slowly-rotating ge¬ 
ometry to second order in the spin. We applied this tech¬ 
nique to study static, axisymmetric tidal perturbations 
of a spinning object. Remarkably, the perturbation equa¬ 
tions can be solved analytically in vacuum. We provided 
the explicit form of the metric describing the exterior ge¬ 
ometry of a spinning object distorted by an axisymmetric 
tidal field to second order in the spin. 

Because of spin couplings, an external quadrupolar 
electric tidal field deforms the quadrupole moment M 2 of 
the central object up to quadratic order in the spin, the 
dipole and octupole current moments Si and S '3 to lin¬ 
ear order in the spin, and also deforms the monopole and 
hexadecapole mass moments Mq and M 4 to quadratic or¬ 
der in the spin. Correspondingly, for a spinning object 
a new class of different Love numbers emerges, while the 
standard Love numbers acquire spin-induced corrections. 

When the central object is a spinning BH, the met¬ 
ric simplifies considerably. Similarly to the nonrotating 
case, we have shown that the multipole moments of a 
Kerr BH are not affected by the tidal field at least up 
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to quadratic order in the spin, and thus the correspond¬ 
ing Love numbers are zero. This implies that the no-hair 
relation (15811 is not affected by perturbative spin-tidal in¬ 
teractions. This result is strongly based on the discrimi¬ 
nation between the solution describing the tidal field and 
that describing the linear response of the central object, 
which we have discussed in some detail. Furthermore, we 
provided the metric describing a tidally-deformed Kerr 
BH in concise form to second order in the spin and com¬ 
puted various physical properties of the solution, includ¬ 
ing the epicyclic frequencies and the intrinsic geometry 
of the event horizon. 

These findings have potentially important impli¬ 
cations for gravitational-wave phenomenology with 
ground-based detectors [^-Q- For example, current 
gravitational-wave templates for compact binary inspi¬ 
rals adopt Love numbers which are valid for nonspinning 
objects, neglecting the spin corrections introduced here. 
In this context, it is crucial to understand how the rota¬ 
tional Love numbers defined here enter the gravitational 
waveforms, similarly to what has been done in the past 
for the nonrotating case [l|, fiol - l^ . 

Our deformed BH solution might be of interest for null 
tests of the Kerr geometry and tests of general relativity 
based on various observations, e.g. tra cing of BH shad¬ 
ows with the Event Horizon Telescope [^, detection of 
quasiperiodic oscillations in the signal emitted by accret¬ 
ing BHs with the X-ray telescope LOFT and ob¬ 

servations of gravitational-wave from extreme-mass ratio 
inspirals with eLIS A [U HO, HH . A more detailed char¬ 
acterization of the solution and its phenomenological ap¬ 
plications are left for future work. 

In a subsequent work [95|, we plan to use the results 
presented here to compute the Love numbers of a spin¬ 
ning neutron star explicitly. This requires solving the 
perturbation equations in the interior of the neutron star, 
taking into account also fluid perturbations and matching 
the interior solution with the exterior metric presented in 
this paper. Another natural question we wish to answer 
is whether the spin-induced corrections to the tidal Love 
numbers of a neutron star satisfy some nearly-universal 
relations as their static counterpart [2^, [s^. 

In this work, we focused on the most relevant, polar-led 
L = 2 perturbations, i.e. we assumed that the tidal field 
has only a quadrupolar electric component at the leading 
order. Extending our results to the axial-led sector (|8]) 
(thus including a magnetic component of the tidal field) 
and to other values of £ is straightforward. As explained 
in the end of Sec. IIV A[ novel families of Love numbers 
will emerge and they would correct the multipolar struc¬ 
ture of a spinning neutron star or of a spinning BH in 
agreement with the selection rules discussed in Sec. HIl 

Indeed, the corrections to the multipole moments of 
a tidally-distorted spinning neutron star will modify the 
approximate three-hair relations that exist for isolated 
compact stars (dsl - l^ . For example, the mass quadrupole 
moment M 2 would acquire 0 (x) corrections proportional 
to the £ = i magnetic component of the external tidal 


field and would acquire O(x^) corrections proportional to 
the £ = 4: electric component of the external tidal field. 
Higher multipole moments would also be modified ac¬ 
cordingly to the selection rules discussed in Sec. HIl Such 
corrections will be explicitly presented in a forthcoming 
publication [^ . 

Although our analysis extends part of the results of 
Ref. to quadratic order in the spin, we did not 
attempt to identify the tidal moments of the source. 
It would be interesting to complement our analysis by 
matching the exterior solution presented here (which is 
valid in the region i? < r <C ro) to a post-Newtonian 
metric describing the source of the tidal field in an over¬ 
lapping region |55| . This would allow us to express the 
constant a in terms of the source parameters up to a 
certain post-Newtonian order. 

A major limitation of our results is the assumption of 
axisymmetry, i.e. we restricted to to = 0 tidal perturba¬ 
tions. The tidal field produced by an orbiting companion 
is not axisymmetric, so our results describe only the ax- 
isymmetric parl[I3 of such tidal interaction (the axisym¬ 
metry of the background guarantees that modes with dif¬ 
ferent TO are decoupled from each other). To second order 
in the spin, a nonaxisymmetric tidal field would induce 
precession of the central object’s angular momentum, in¬ 
troducing time dependence in the problem. This effect 
can be avoided by restricting the analysis to first order 
in the spin, as done in Refs. [H, In this case our 
method can be straightforwardly extended to compute 
the new Love numbers in the to ^ 0 case to linear order 
in the spin. In principle, however, our method is also well 
suited to deal with time-dependent perturbations and we 
hope to come back to this interesting generalization in the 
future. Indeed, the slow-rotation framework presented in 
this paper can be applied to the case of time-dependent 
perturbations with only minor modifications [72|. Such 
extension can be also relevant for a variety of important 
problems, including the computations of the quasinormal 
modes of spinning neutron stars [9^ to second order 
in the spin. 

Note added. After this work was in its final stages, 
we were informed that Poisson and Landry are extend¬ 
ing the results of Ref. to the case of a material 
object and of nonaxisymmetric tidal perturbations. Our 
formalism and theirs are in various ways complementary 
to each other: differently from Ref. [63 we extended the 
perturbative analysis to second order in the spin, but we 
restricted to the axisymmetric case and did not include 
a quadrupolar magnetic component of the tidal field. 


It is important to stress that the Love numbers do not depend 
on the details of the source. Therefore, our results describe the 
(axisymmetric, quadrupolar) rotational Love numbers of a spin¬ 
ning body immersed in a generic tidal field, regardless of the 
symmetry of the latter. 
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Appendix A: Background geometry to second order in the spin 

Here we give the metric coefficients entering the background geometry ([2|) to second order in the spin [tI, [t^ : 




mo(r) = ( Sm — 


m3 


m2{r) = -x^ 


]VP 


SmM 


r-2M r3(r-2M )J ’ 


77T,2(r) =-;;—(5M — t) (r — 2M) 


2Mh 

j 2 (r) = 


+ 


M(M - r) (3r2 - 2M^ - 6Mr) -f 3r^(r - 2Mf tanh“ 


M^X^iM + r) 


M 


r — M 


2M2r(r - 2M) 
V 2 ir) = - 


^M{r — M) (3r^ — 2M^ — 6Mr) — 3r^(r — 2M)^ tanh ^ ^ 
r4 M ^ {r-MJ r{r-2M) j ’ 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


the meaning of the various constants is explained in the main text. Note that we have factored out the spin dependence 
of the second-order terms, so that 5m and 5q are 0{x^) numbers which multiply terms proportional to x^- 


Appendix B: Stationary perturbations of a slowly-rotating relativistic object 


In this Appendix we briefly present the derivation of the field equations; to reduce the risk of typographical errors and 
facilitate comparison with our results, we made the entire calculation available online in the Supplemental Material. 

As a background, we consider the spinning geometry Q to second order in the rotation rate (cf. Appendix]^ and 
we consider a harmonic decomposition of the metric perturbations as 


A r, <fi) + r, z?, (/?) (Bl) 

with 

Ks$ \ 

hisi h[Sf^ 

hi sin 
* hi sin 


SgT. = 


/o 0 
* 0 
* * 
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(B2) 
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* airPiY^ 

dpi 

PYP 


(B3) 

* * 

r2 + G^W^] 

r^G^X^ 



* 

r^ sin^ -d [X^y^ - 

G^IU^jy 



where asterisks represent symmetric components, 1 /g^r — 1 — 2A4 jr, = Y^{'d, ip) are the scalar spherical 

harmonics and we have defined 




(X^ W^) = 1 - cot - cot ^y^^ - ^ 1 , 


y^ 

r£ ,VV 


(B4) 

(B5) 


which are related to the vector and tensor spherical harmonics, respectively. Here and in the following, a sum over the 
harmonic indices i and m (such that |to| < £) is imolicilF^. Under parity transformations {'& —> tt— z?, p —>■ p+n): polar 
and axial perturbations are multiplied by (—1)^ and (—respectively. The odd and even sectors are also referred 
to as “axial” and “polar” and we shall use the two notations indistinctly. The functions {Ho,Hi,H 2 , K, G, r]o,ViY and 
{ho, hi, h 2 Y only depend on t and r and describe the polar parity metric perturbations and the axial parity metric 
perturbations, respectively. In the following, we adopt the Regge-Wheeler gauge [1^ and set gf = = h^ = 0 

through a gauge choice. 

Using this decomposition, we can solve vacuum Einstein’s equations perturbatively in the spin. Because of 
the transformation properties of the perturbation functions, the linearized equations naturally separate into three 
groups By denoting the linearized Einstein equations as = 0, the first group reads 


SS(^i) = {A[^^ + cosd + cos^ 'd)y^ + + B^P cos'd) sin-dU!^ = 0, 


(B6) 


where I = 0,1, 2, 3 corresponds to 6 £u = 0, 5£tr = 0, SSrr = 0 and S£m + 6 £,p,p/sm'd^ = 0, respectively. The second 
group reads 


yf 

5£(^l^)= {aP + aP cosz? + aP cos^ - {PP + pP cosi? + Pp cos^ 

+ (gP + fjP cos'd) sini?y^ + {^P + cos'd)X^ + xP sin^JIU^ = 0, (B7) 

Y^ 

{pP + pP cos'd + pP cos^ 'd + aP sin^ cos'd + cos^ 'd + A^ sin^ -d) . 

+(d^^ + cos 'd) sin ^Y^ + xP^X^ - + PP'’ cos i?) sin dW^ = 0, (B8) 

where L = 0,1 and the first equation corresponds to S£t^ = 0 and d£r^ = 0, whereas the last equation corresponds to 
S£tip = 0 and S£r(p = 0. Finally the third group is 




= {fi + cos'd) sin-dy;^ + {gt, + gt cosd)y^ + ki sin^ 'dY^ + {sc + S£ cos^ -d )-—- + {tc + tc cos^ 'd)W'^ = 0 


sini? 


(B9) 


(!£■(_) = {gc + gt cos'd) sin'dy;^ - {fc + fi cos'd)Y‘c + kc sin^ 'dY'^ - {tc + U cos"' -d )-—- + (s^ + se cos^ 'd)W^ = 0, 

’ sinv 


(BIO) 


corresponding to 6 £^,p = 0 and 6 £^^ — 6 £,p,p/sm'd^ = 0, respectively. In the equations above, X^ and are the 
tensor spherical harmonics defined as in Eq. (IB5I) . The coefficients appearing in these equations are Imear and purely 
radial functions of the perturbation variables. Furthermore, they naturally divide into two sets accordingly to their 
parity: 
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Axial: 
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Furthermore, from now on we will append the relevant multipolar 
index i to any perturbation variable but we will omit the index m, 
because in an axisymmetric background it is possible to decouple 


the perturbation equations so that all quantities have the same 
value of m. 
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The explicit form of the coefficients is given in the online notebook in the Supplemental Material. Note that the 
coefficients above are purely radial functions, i.e., the entire angular dependence has been completely factored out in 
the linearized Einstein equations. 


a. Separation of the angular dependence 


The decoupling of the angular dependence of the Einstein equations for a slowly-rotating star to first order in the 
spin was performed by Kojima [sB] (see also [65|) by using the orthogonality properties of the spherical harmonics. 
The procedure has been extended to the case of slowly-rotating BHs to second order in Refs. 0, [^. Here we adopt 
the same technique. 

The decoupling is achieved by computing the following integrals 


0 = y 56(^1) , 

’ sin 17 

0 = / dn- - \ - 

0 = / dfl- - \ - 

J £(£+!)-2 


^ ' Sin 17 ^ ^ 


^ ' Sin 17 ^ ^ 



(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 


for I = 0,1, 2,3 and L = 0,1. These integrals might be evaluated analytically for generic {i, m) by using the properties 
of the spherical harmonics [H, [t^, [zll • The resulting equations have the form dl])-®. For simplicity, here we perform 
the integrals above explicitly. For example, for i = m = 2 one obtains the system of equations 
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(B16) 

(B17) 

(B18) 

(B19) 

(B20) 


which describe the full set of £ = m = 2 perturbations to second order in the spin. Similar equations can be found 
for other values of £ and \m\ < 1 . 

To summarize, our decoupling procedure in the slow-rotation limit allows us to obtain a system of 10 coupled, 
ordinary differential equations. The mixing of different angular functions in Eqs. (IB6|) - (IB10I) yields a mixing of 
perturbation functions with multipolar indices l-\-\ and £ — 1 to first order in the spin and with multipolar indices 
i -\-2 and £ — 2 to second order, respectively. The explicit form of the final radial equations is available online in 
the Supplemental Material. 
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1. Perturbations with 1= \ and I = Q 


In the derivation presented in this Appendix we have always assumed (. > 2. However, in the rotating case axial 
and polar perturbations with i = 2 are respectively coupled to polar and axial perturbations with I = ?> and £ = 1 
and with axial and polar perturbations with £ = 0 and £ = 4. Perturbations with £ = 0,1 satisfy a different set of 
equations than the one presented above. The perturbation equations for £ = 1 can be found by neglecting Eqs. (IB14I) 
and (IB15I) (which are clearly ill-defined when £ = 1) and by using a residual gauge freedom in order to set = 0 
and = 0 in the ansatz for the metric llOCll |. For £ = 0, only scalar polar perturbations exist and one can set 
^'(0) = Q = Q without loss of generality. 


Appendix C: Sources for the inhomogeneous equations governing tidal deformations 


The source terms entering the axial perturbations for L = 3 and L = 1 read 





(Cl) 

(C2) 


whereas the source term entering Eq. dni) reads 
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The source term entering Eq. m reads 
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Finally, the source terms entering Eqs. dsnii-dii]) read 
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Appendix D: Explicit solution for a tidally-deformed spinning vacuum geometry to first order in the spin 


Here we present the explicit solution for a tidally-deformed spinning vacuum geometry to first order in the spin. To 
zeroth order in the spin and in the tidal field, the background solution was given in Appendix]^ For L = 2 polar-led 
tidal perturbations, the nonvanishing metric components to zero order in the spin and to first order in the tidal field 
are given in Eqs. (fT^-(fn|. Finally, to first order in the spin and in the tidal field, the only nonvanishing metric 
functions are the following L = 3 and i = 1 axial components: 
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where again y = r/M and 01^3 and 71^3 are integration constants related to the homogeneous problem associated 
with Eqs. (HU) and (1171) (cf. main text, Eq. (1281) and Table U). An analytical solution can be found in closed form 
also to second order in the spin, by solving Eqs. dni), (HU), dini) and (I 2 T]) analytically, but it is too cumbersome to be 
displayed here. The explicit form is given in a notebook in the Supplemental Material. 

The hrst-order solution generalizes that found by Poisson (55j | to the case of generic vacuum and reduces to it in 
the BH case with m = 0 and with a purely electric quadrupolar tidal field. To linear order in the spin, the coordinate 
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transformation that brings our BH metric to the axisymmetric, electric-led solution found in 
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where dr/dr* = 1 — 2M/r and the term proportional to the light-cone time v in tp is needed to eliminate the gauge 
term proportional to in Ref. [s^. In the notation of Ref. [s^, £q = M and x are the same quantities as 

the ones we used, and we remark that the normalization of the spherical harmonics in Ref. differs from ours. 


14 Actually, our solution coincides with that found in [H after Poisson and Phil Landry for clarifying this point in a 

adding a term ^ in the function defined in Table III of private communication. 

Ref. Such term, however, is arbitrary as it can be reab¬ 

sorbed in a shift of the body’s angular-momentum vector. We 
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Appendix E: Low-velocity expansion of the energy change AE 

For our tidally-distorted vacuum solution, the low-velocity expansion of the energy change AE per logarithmic 
interval of the orbital frequency, Eq. m, reads 
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By comparing Eq. (lEll) with Eq. dm, we can identify the first multipole moments, Eqs. (l44D -(l48 p . Although not dis¬ 
played here, we have also computed the coefficients of and and checked that they agree with the expansion (1411) 
once the dehnitions (H^ -(l48 p are imposed. This is a consistency check of our solution. 
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Appendix F: Explicit solution for a tidally-deformed spinning BH to second order in the spin 


When the central object is a BH, regularity at the horizon simplifies the metric perturbations considerably and the 
full solution can be presented in compact form. The nonvanishing components of the line element ([2]) in the case of 
a spinning BH deformed by a stationary, axisymmetric quadrupole-led tidal field to second order in the spin read 
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where y = r/M. This solution generalizes that found by Poisson to quadratic order in the spin, but restricting to 
the axisymmetric (m = 0) case with a purely electric quadrupolar tidal field. The coordinate transformation between 
our solution at that found in Ref. [s^ is given in Appendix 
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